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ABSTRACT 

We study the collapse and virialization of an isolated spherical cloud of self-gravitating 
particles initially at rest and characterised by a power-law density profile, with expo- 
nent < a < 3, or by a Plummer, an Hernquist, a NFW, a Gaussian profile. We find 
that in all cases the virialized structure formed after the collapse has a density profile 



decaying, at large enough radii, as 



and a radial velocity dispersion profile de- 



caying as ~ r . We show that these profiles originate from the physical mechanism 
responsible of the ejection of a fraction of cloud's mass and energy during the collapse 
and that this same mechanism washes out the dependence on the initial conditions. 
When a large enough initial velocity dispersion is given to the cloud particles, ejec- 
tion does not occur anymore and consequently the virialized halo density and velocity 
profiles display features which reflect the initial conditions. 

Key words: Galaxy: halo; Galaxy: formation; globular clusters: general; (cosmology:) 
dark matter; (cosmology:) large-scale structure of Universe; galaxies: formation 
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1 INTRODUCTION 

It is well known that some of the properties of virial- 
ized structures formed through the gravitational collapse 
of an isolated system depend on the initial conditions 
in various ways; however, some other statistical fea- 
tures of these structures seem to be generated by the 
dynamical mechanis ms acting d uring the relaxation 
proces s (see 
19881); iTheis 



(2006): 



g., Ivan Albadal (ll982Tl; lAarseth et ail 



Spur zed (1 19991): iBoilv 



Joyce. Marcos fc Svlos Labini 



Athanassoulal 
( 2009al lbl)~ 



IVisbal. Loeb fc Hernquistl (|2012l ) and references therein) 



For instance, when a spatially uniform and isolated 
cloud of density po and with a small enough initial ve- 
locity dispersion collapses, it forms, in a time-scale td = 
\J 3n I (32G po) , a virialized structure characterised, for radii 
larger than the size of a dense inner core r c , by (i) a power- 
law density profile n(r) ~ r~ 4 , (ii) a radial velocity disper- 
sion which decays in a Keplerian way a^.{r) ~ r _1 and (iii) 
a pseu do phase-space de nsity that decays as n(r)/o~r(r) ~ 
r" 5/2 (|Svlos Labinill2012l) . 

These behaviors were shown to be associated with the 
ejection of a fraction of the initial mass and energy from 
the system during the collapse: in particular, particles with 
binding energies close to zero follows almost radial orbits 
around the dense core, thus displaying a Keplerian veloc- 
ity dispersion. Considering that they conserve their total 
energy it is possible to show that the density profile fol- 



lows the observed n(r) ~ r~ 4 decay (|Svlos Labinil [20121 ) . 
Note that the crucial role of particles ejection for the for- 
mation of non-trivial density and velocity profiles has not 
been p roperly appreciated in the l i teratur e (see discus- 
sion in Ijovce. Marcos fc Svlos Labinil (|2009al ) ; ISvlos Labinil 
(EoH)). 

In order to study whether the behaviors (i)-(iii) above 
are found in the collapse of spherical clouds with dif- 
ferent initial density profiles than uniform and to deter- 
mine whether in more general cases there remains mem- 
ory of the initial conditions after the collapse phase, we 
consider (i) a set of clouds with an initial density profile 
of the type n(r) — n a /r a , where n a is a constant and 
^ a < 3 and (ii) density profiles without a simple scale 
free behaviour, e .g., a Plummer sphere, an Hernquist profile 
(|Hernauistll990l). a Navarro Fre nk and White (NFW) profile 
(|Navarro. Frenk fc W hite 1997) and a Gaussian profile. 

Note that while the formation and evolution of self- 
gravitating dissipation-less systems is a theoretical interest- 
ing and complex problem, the simple models considered here 
may only help to develop tools and analysis to understand 
more realistic systems for which it is then necessary to con- 
sider the role of merging, of dissipation and the fact that 
they cannot be simply considered isolated. Whether these 
complications qualitatively change the evolution with re- 
spect to the simple systems considered here will be stud- 
ied in a forthcoming publication. However, the collapse of 
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an isolated cloud of self-gravitatin g particles may describe 
the formation of elliptica l galaxies (|Henonlll964l ; lvan Albadal 
ll982l : ISvlos Labinill2012l ) which are characterised by the de 
Vaucouleurs R 1 ^* law. A possible analytical approximation 
to the de-projec t ed R 1 ^ 4 is given by the Hernqu ist profile 
|Hernauisdll990l ; IVisbal, Loeb fc Hernouistl 120121 ) which is 
indeed characterised by a r -4 decay at large radii. 

The paper is organised as follows. We present in Sect [2] 
the details of the numerical simulations that we have per- 
formed. The case of an initially scale free profile is discussed 
in Sect(3] We then consider in Sect[4]more complex density 
profiles, which have a characteristic length. Finally in SectfS] 
we discuss the results draw our main conclusions. 



2 THE SIMULATIONS 

The initial conditions are generated by distributing N par- 
ticles randomly with density n(r), where r is the distance 
from the origin of coordinates. For the scale free density pro- 
files n(r) = n a r~ a , where ^ a < 3: for larger values of 
a and r — > there are too many nearest neighbours with 
arbitrarily small distance. In such a situation an accurate 
numerical integration of the equations of motion becomes 
very problematic. We put a lower and an upper cut-off: the 
latter is of the order of the softening length of the code and 
the former, Ro , is larger than any characteristic scale of the 
profile (if present). 

We generate a series of spherical clouds of particles, 
with N — 10 4 and with different initial virial ratio bo = 
b(t = 0) = 2K /W a (where K/W is the kinetic/potential 
energy). In order to assign velocities we take the velocity 
components to have a uniform probability density function 
(PDF) in the range [— Vb, Vo] so that the modulus of the 
velocity is constrained to be in a sphere of radius Vo. The 
velocity PDF is thus 

3 



9(v) 



for v ^ Vb 



(1) 



and zero otherwise. The mean square value of the velocity 
is 

/•Vfa 

2\ 



2 / u 3 2 Wobo 
vg(v)dv=-V = 1 ^. (2) 



We used the publicly availab le tree-code GADGET 
(jSpringel et al.ll200ll : ISpringell[2005l ) to run the simulations. 
We have c hosen the same com bination of numerical param- 
eters as in ISvlos Labinil (|2012f l, to which we refer the reader 
for further details on the numerical issues. 



3 SCALE-FREE PROFILES 

We first discuss the relaxation of a spherical and isolated 
cloud, with a with a power-law profile, initially at rest. Then 
we consider the effect of a non-zero initial velocity disper- 
sion. 



3.1 Virial ratio 

The virial ratio for all system particles b a (t) and the virial 
ratio for bound b n (t) have a different amplitude (see upper 
left and right panels of Figfl}: this can be easily explained as 




Figure 1. Upper left panel: Virial ratio of all particles. Upper 
right panel: Virial ratio of particles with negative total energy. 
Bottom left panel: Fraction of particles with positive energy. Bot- 
tom right panel: Behaviors of the (normalised) energy exchanged 
Ajv(i) (see text for details). 



due to the fact that part of the particles have been ejected 
during the collapse. Indeed, the fraction of particles f p {t) 
with positive energy has an abrupt change becoming larger 
than zero for t 55 rjj (see FigfTJ — bottom left panel). For 
a ^ 2 only f p (t) w 1% of the particles gain positive energy 
during the collapse while for a = we find f p (t) > 20%: 
the amount of ejected particles decreases with a and so the 
difference between btot(t) and b neg {t). 



3.2 Ejection 

As dis c ussed in detail in Ijovce. Marcos fc Svlos Labinil 
|2009al ); ISvlos Labinil (|2012i ) the mechanism allowing par- 
ticles to gain energy during the collapse phase arises from 
the interplay between the growth of density perturbations 
in the collapsing cloud and the finite size of the system. In 
practice, particles that were originally in the outer shells of 
the cloud arrive later than the others in the centre, thus mov- 
ing, for a short time interval around ru, in a rapidly varying 
gravitational potential. As noticed above, even when a > 
some particles acquire positive energy at about t » to, i.e. 
when all other particles have already arrived to the centre 
and inverted their velocities. The difference between the case 
a = and a > 0, and so the variation of f p with a, is due to 
the different time of arrival at the centre of particles initially 
placed at different radii. 

In the fluid limit and neglecting density perturbations, 
a particle initially placed at ro obeys to the following para- 
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metric equations: 

r(0 = y(l + coB(fl) 



(3) 



*«) = 
where 

M(ro) = 



8GM(r 



(f + Bin(0) , 



mn a r a 4irr dr ■ 



47T 

3-a 



mn a r 



(4) 



Thus we get that a particle initially at ro arrives at the origin 
for £ = 7r and thus at 



r(r ; a) = 



8GM(r ) 



(3 — a)7rro 
32mn a G 



(5) 



i.e., particles initially placed in the outer shells arrive later 
than the others for a > even if they follow their unper- 
turbed trajectories. Instead, in this same limit, for a — 
all particles should arrive simultaneously at the centre, i.e., 
r(ro, 0) = td Vr . 

Noticing that for ro = -Ro we have M(Ro) = miV and 
po = M(R )/V(R ) = 3miV/(47r7?g) we find that 



r(r < -Ro; a) < r(R ; a) = t d Va 



(6) 



Thus the maximum time taken by a particle to arrive to the 
centre does not depend on a. Given that we detect particles 
with positive energy only for t ^ td, being td the time taken 
by the particles in the outer shell to arrive at the centre, 
we may conjecture that, as it occurs for the a = case, 
ejected particles are those which initially lie in the outer 
shells: indeed Figd] (upper right panel) shows that this is 
the case. 

In order to estimate the scaling of the fraction of ejected 
particles as a function of a, let us suppose that all ejected 
particles belong initially to the outer shell, i.e. they have 
at time t = radius ro 6 [Ro — AR, Ro]. The number of 
particles in this last shell scale as Q 



AiV : 



Aim a r 2 a dr = N I 1 



ARY 
Ro J 



where we used the normalisation condition (N is fixed) 
_ (3 - a)N 



(7) 



(8) 



The scaling behaviour of the fraction of ejected particles as a 
function of a derived from Eq[7] i.e. f p = AN/N, is shown in 
Fig[2] (upper left panel), together with the simulation points. 
(Note that the amplitude is fixed by the choice of AR.) One 
may see that a very rough agreement is found: in principle, 
one should use the true probability p(ro) for a particle to 
be ejected as a function of its initial radial position ro- For 
in stance, by assuming p(rp) ~ r„ with /3 > 3, as measured 
in ljovce. Marcos fc Svlos Labinil (|2009aT ) for the case a = 0, 
we find 



AN 

IT 



3-a + /3 



(9) 



This corresponds to assuming a radial probability distribution 
for a particle initially placed at ro for being ejected at t > td of 
the type p(ro) = ©(ro — -Ro + AR) where 0(x) is the Heaviside 
step function. 



r J R o 



0.5 1 1.5 2 2.5 3 0.2 0.4 0.6 0.8 1 




Figure 2. Upper left panel: Number of ejected particles as a 
function of their initial radial position for a = 0.5 and a = 2. 
B Upper right panel: Fraction of the particles with positive en- 
ergy as a function of a. The dashed line corresponds to a radial 
probability of ejection described by the Heaviside step function 
(i.e., p(ro) = ©(ro — -Ro + AR) — see EqO while the solid line 
to a power-law radial probability p(ro) ~ with /3 = 4 (i.e. 
Eq|9]l. Bottom panel: Behaviour of the energy e p (t) for some of 
the ejected particles as a function of time, normalised to the en- 
ergy at t = 5te> and for a = 0. 



Even in this case the agreement is very rough: however the 
assumption that p(ro) has the same exponent for all values 
of a is probably too strong. Indeed, as mentioned above, the 
physical mechanism of ejection is due to a complex interplay 
between the fluctuations present in the distribution and the 
finite size of the system. 



3.3 Energy exchange 

It is interesting to consider how individual particles gain en- 
ergy as they pass through the core. In Fig[5] (bottom panel) 
it is plotted the behaviour of the energy e p (t) for some of the 
particles that will be ejected after the collapse as a function 
of time, normalised to the energy at late times (i.e. t = 5td) 
and in the case a — 0. One may note that there is a single 
abrupt change of the energy at t ~ td, a fact that shows 
that particles are being accelerated by large scale changes of 
the mean field potential, instead of getting rapid kicks from 
encounters with individual particles. 

A statistical manner to monitor how particles energy 
changes during tim e consists in measuring the quantity 
l|Svlos Labinill2012t ) 



<A 2 (t)> 



(e* (i) - e p {t)f 



N{N - 1) 



<e(i)> 2 



(10) 
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where ei,(i) the average energy per particle is defined as 

e£\4(*) 



(e(t)> 



N 



(11) 



Fig[T] (bottom right panel) shows the behaviour of A%(t) = 
(A 2 (t)) - (A 2 (0)) in the various simulations: this changes 
with time when particles exchange energy, i.e. when the par- 
ticle normalised energy distribution changes. 

It is interesting to note that while for ot = the ex- 
change of energy occurs simultaneously with the strongest 
phase of the collapse at t ~ td , for a > particles exchange 
energy also for t < td. This occurs because, for a > the 
i th particle, initially located at r l > rjj S> 0, arrives at the 
centre, if it follows its unperturbed trajectory, later than 
the j th one (see Eq(SJ. Thus when the i th particle arrives 
at the centre, inverting its velocity, it moves in a varying 
potential field generated by all other particles which were 
initially placed at distance < rj. However the exchange of 
energy, as long as ro < Ro, is not enough to allow a frac- 
tion of the particles to escape from the system. This is be- 
cause the gain of kinetic energy depends on the deepness of 
the potential well and thus on the to tal mass generating it 
(| Joyce. Marcos fc Svlos Labinill2009al '). 

Finally, we have plotted in Fig[3]the average radial ve- 
locity at time t of the particles that at later times, i.e. at 
5td, have positive (+) and negative energy (-). One may 
notice that the change in the radial velocity direction oc- 
curs sooner for the latter particles than for the former ones. 
Particles that will be ejected arrive later than all others, 
having about the same velocity independently on a. Indeed, 
as they belong to the outer shell, they are attracted by the 
same gravitational field. For this reason the dynamics of the 
ejection is the same for different values of a: what changes 
is, as previously discussed, only the amount of mass and 
energy that is ejected from the system. 



3.4 Density profile 

In the case of an initially uniform density profile, i.e. a = 0, 
the radial density profile of the virialized structure formed 
by bou nd particles after the collapse was found to have a r~ 4 
decay . |jovce. Marcos fc Svlos Labinil l2009al ; ISvlos Labinil 
I2012T I. For a > one may notice (see FigfJ] — upper left 
panel) that at small radii the different profiles have differ- 
ent slopes which are close to the initial ones (for a value of 
the radius of order of the force softening length e the profile 
seems to flatten, probably due to the fact that the interac- 
tion is smoothed) A simple interpretation of this behaviour 
is the following: particles in the dense core were already very 
strongly bound at t = 0, as the density has a cusp for a > 
by construction; thus the energy exchange during the violent 
phase of the collapse represents a small perturbation that is 
not able to alter significantly the inner particles trajectories. 

On the other hand for r/Ro ~ 1 all profiles converge to 
the same shape, i.e. n(r) ~ r ,_4 Va. We thus conclude that 
the dynamical mechanism underlying the formation of the 
r -4 should be the same for all a values, i.e. it is related to 
process giving rise to mass and energy ejection. 

A simple way to focus on the properties of the r -4 tail is 
to monitor how this changes during time evolution. Indeed, 
as mentioned above, in the case a = 0, the large scale tail 
which decays as n(r) ~ r~ 4 is made of particles with almost 




Figure 3. Behaviour of the average radial velocity for particles 
that at t = 5t£) are bound (-) and ejected (+). 
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Figure 4. Upper left panel: Density profiles n(r) at t = 5t£>: the 
amplitude for a > has been arbitrarily normalised to have the 
same large radii amplitude of n(r) for a = 0. Upper right panel: 
Square value of the radial component of the velocity (v^(r)) as 
a function of the radius at t = 5td for different values of a. 
Bottom left panel: Pseudo phase-space density $(r) as a function 
of the radius at t = 5t£> for different values of a: the different 
curves have been normalised to have same amplitude at large 
radii as the a = case. Bottom right panel: Average energy per 
particle (e p (r)) (absolute value) as a function of the radius and 
for different values of a. 
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zero energy orbiting around the dense core at smaller radii. 
As time passes these particles may reach larger and larger 
radii and thus the tail extends to larger and larger radii. This 
is indeed what is observed not only in the case a — but also 
for a > 0. This evidence corroborate the conclusion that the 
density profile tail properties have the same physical origin 
for all a values. We will come back on this point in what 
follows. 



3.5 Radial velocity dispersion profile 

In FigfJ] (upper right panel) it is shown the mean square 
value of the radial velocity as a function of the radius 
cr 2 (r) = (Vr(r)} at t — 5td for different values of a. One 
may note that the different curves are almost indistinguish- 
able and that at large radii, i.e. r/Ro > 0.1, they show an 
approximated Keplerian <r 2 (r) ~ r _1 behaviour. 

The time evolution of (v^(r)) is again very similar for 
different valu es of a and compatib le with the physical model 
presented in ISvlos Labinil (|2012h : a few particles with al- 
most zero energy continue to increase their distance from 
the structure core, remaining however bound. Note that, at 
a given time, the break of the r~ l behaviour is found at 
the same radii for different values of the a: this reflects the 
fact that particles belonging to the tail not only obey to 
the same equation of motions but also are subjected to the 
gravitational effect of approximately the same core mass. 



3.6 Pseudo phase-space density 

The behaviour pseudo phase-space density $(r) = 
n(r)/o~r{r) (see Fig(4] — bottom left panel) is a combina- 
tion of the density and radial velocity dispersion profiles. In 
agreement with the previous results, at large enough radii it 
shows the same behaviour independently on a; indeed, a fit 
with a function of the type <J?(r) ~ — ^- ~ r _5//2 describes 
only the large scale tail. Finally we note that, as for the 
radial velocity dispersion, at a given time, the break of the 
r~ 6 / 2 tail is found at the same scale for different values of a 
(in this, at t = 5td, there is a bump at r/Ro « 3). 



3.7 Energy distribution 

The distribution of particles total energy at t « 5td, i.e. 
when the virialized structure is already formed, is shown in 
FigE] (left panel). One may note that the larger is a the 
smaller is the fraction of particles with positive energy, a 
result clearly in agreement with the ejection mechanism dis- 
cussed in Sect 13. 2 1 

The average energy of a particle with radial distance 
r from the structure centre is displayed in Fig(4] (bottom 
right panel). At large radii, i.e. r > 1/10-Ro, the energy 
displays a r ,_1 behaviour, in agreement with the fact that the 
gravitational potential energy is dominated by the inner core 
and thus that the velocity dispersion presents a Keplerian 
decay. On the other hand, at small radii, the different density 
profile give rise to slightly different potential and kinetic 
energy profiles. 
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Figure 5. Left panel: Particle total energy distribution for dif- 
ferent values of a Right Panel: Velocity anisotropy at t = 5td- 
When /3(r) > the velocity dispersion is dominated by its radial 
component. 



3.8 Velocity anisotropy 

FigE] shows the behaviour, in the different simulations, of 
the velocity anisotropy parameter defined as 



/3(r) = 2 - 



vt(r) 
v r (r) 



(12) 



where vt(r)/v r (t) is the velocity in the transversal/radial di- 
rection. In the outer regions of the virialized structure where 
the density scales as n(r) ~ r~ 4 , rad ial orbits dominat e 
the motion, as it was already noticed by van Albadal(ll982tl : 
iTrenti. Bertin fc van Albadal (|2005l ): ISvlos Labinil (|2012h for 



the a = case. It is interesting to note that not only the 
behaviour of f3{r) is essentially the same for different values 
of a, but also its amplitude. This, together with the evi- 
dences already discussed from the radial velocity dispersion 
and its time dependence, confirms the fact that the particle 
dynamics in the large scale tail is the same for all values of 
a considered. 



3.9 The case of a nonzero initial velocity 
dispersion 

For an initially uniform cloud (i.e., a — 0) and for warm 
enough initial velocity dispersion e jection of ma s s and en- 
ergy does not occur anymore. In ISvlos Labinil (|2012T ) we 
found that the critical value of the initial virial ratio was 
ba ~ —1/2: for bo > &o ejection occurs and for bo < bo it does 
not occur anymore. The value of bp is approximate as it was 
conjectured to depend on the number of points N used in 
the simulation. The reason for not having ejection anymore 
is due to the fact that an high enough velocity dispersion 
cause shell crossing well before the collapse time td- Such a 
shell crossing erase the temporal lag of the outer shell, thus 
avoiding that particles initially placed at the outer bound- 
aries of the system arrive later, then moving in a rapidly 
time varying gravitational field. 

A similar mechanism should be working also for a > 
0. For each value of a and for each TV one should be able 
to find a critical value of the initial virial ratio, i.e. b§ = 
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Figure 6. Density profile and velocity dispersion for a = 1 (upper 
panels) and a = 2 (bottom panels) and bo = 0, —1/2, — 1 at t = 
5td. 

&o(-/V; a). Here we have limited our study to bo = —1/2, — 1 
and a = 1, 2 for which we have not found mass and energy 
ejection. Instead of determining the precise value of &§ we 
are interested here in the shape of the density and velocity 
profiles when ejection does not occur. 

In this case the collapse is more gentle, and the virial ra- 
tio shows a series of damped oscillations which drive the sys- 
tem toward the quasi-equilibrium state. In Fig[S]we show re- 
spectively the density profile, velocity dispersion and pseudo 
phase-space density for a = 1, 2 and bo — 0, —1/2, —1. One 
may note that (i) the smaller is bo and the smaller the differ- 
ence between the input and output density profile, (ii) For 
different values of 6o there is a difference in the behaviour 
of all density and velocity profiles, a fact that points toward 
an important role of the initial conditions, (iii) The case 
a — 2 presents the smaller change in the density profile n(r) 
between the initial and final condition: the larger is a and 
the deeper the potential well at small scales and thus the 
smaller the perturbation introduced by a given amount of 
kinetic energy. We thus conclude that when the collapse is 
violent, for the effect on an initial velocity dispersion, there 
is a clear dependence of the virialized structure properties 
on the initial conditions. 



4 RELAXATION OF MORE COMPLEX 
DENSITY PROFILES 

In this section we have considered and initially perfectly cold 
cloud of self-gravitating particles, whose profile is described 
by: (i) the Plummer sphere 

A 
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-■ — ■ Plummer 




■ — ■ Hernquist 




• — » Gaussian 




* — * NFW 
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Figure 7. Upper left panel: Virial ratio of all particles. Upper 
right panel: Fraction of particles with positive energy. Bottom 
left panel: density profiles n(r) at t = 5tq Bottom right panel: 
radial velocity dispersion (v^,(r)) at t = 5to 



(ii) the Hernquist profile (7 = 3), (iii) the NFW profile (7 
2) 

,1 



n(r) 



and the (iv) Gaussian profile 



n(r) = Aexp 





\ 7 


r(lH 









(14) 



(15) 



n(r) 



(' + #)' 



(13) 



where in Eas llSlfT^I A and r s are two constants that are fixed 
by the number of points N p initially belonging to the cloud 
and by the requirement that the radius enclosing half of the 
mass is the same for all profiles, so that the collapse time is 
about the same in the different cases. 

Results are shown in Fig[7] we have normalised the time 
to the characteristic time t 3 for collapse for each density 
profile and the radius to r s defined in Eqs ll3ll5l One may 
note that there are similarities with the behaviours found in 
Figs lH4l In particular in all case it is observed ejection of a 
fraction of the particles and the formation of a n(r) ~ r -4 
density profile and a radial velocity dispersion decaying as 
r . Qualitatively the dynamics is thus similar to the case 
discussed in the previous section. 



5 DISCUSSION AND CONCLUSIONS 

The collapse of an isolated spherical cloud of self-gravitating 
particles represents a paradigmatic example of the dynamics 
giving rise, through a global relaxation process, to a virial- 
ized structure in a quasi-equilibrium state, i.e. a state that 
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has a life-time that diverges with the number of system par- 
ticles and that is far from thermodynamical equilibrium. If 
the collapse is violent enough part of the initial mass and 
energy of the cloud is ejected in the form of free particles. 
This ejection occurs, for the case of an initially uniform 
cloud, when the in itial velocity dispersion is small enough 
(|Svlos Labinill2012h . In this case some non-trivial density 
and velocity profiles have been observed: (i) the radial den- 
sity decays as n(r) ~ r~ 4 at large enough radii, (ii) the radial 
velocity dispersion decays as a 2 ~ r _1 and thus (iii) the 
pseudo phase-spac e density decays as n(r)/af(r) ~ r - -5 / 2 
(jSvlos Labinil2012r ). All these behaviors can be understood 
in the framework of a simple dynamical model in which the 
bound particles that have the largest energies, belong to the 
tail of the profile and are orbiting in a central gravitational 
poten tial generated by the inner dense core (|Svlos Labinil 
l2012h . 

The problem that we have considered in this paper con- 
cerns whether the same density and velocity profiles are ob- 
served when the initial conditions are not represented by a 
uniform cloud of particles. In particular, we have firstly con- 
sidered an initial density profile of the type n{r) ~ r~ a with 
^ a < 3. When the initial velocity dispersion is zero, the 
situation for ot > is similar to the uniform case, i.e. a = 0, 
in various aspects. During the collapse phase a fraction of 
the particles gain enough kinetic energy to escape from the 
system, i.e. there is ejection of mass and energy. The viri- 
alized structure formed after the collapse presents, at large 
enough radii, the same behaviors (i)-(iii) discussed above for 
the initially uniform cloud. In addition, radial orbits domi- 
nate the motion of particles in the density profile tail in all 
cases. The time evolution of these quantities shows that the 
underlying dynamics determining the large radii properties 
of the virialized structures is identical to the a — case. 
These common features are thus generated by the dynam- 
ical mechanism acting during the collapse, and responsible 
of the mass and energy ejection. For this reason they are 
independent on the initial properties of the self-gravitating 
cloud. 

On the other hand, when the initial velocity dispersion 
is high enough, we find that there is no ejection of mass 
and energy so that the density and velocity profiles do not 
present any power-law behaviour at large radii. In this case 
the collapse is more gentle and consists in a few oscillations 
that drive the system to a quasi-equilibrium configuration 
and the virialized structure displays features that depend on 
the initial conditions. 

In addition, we have considered other initial density pro- 
files which are not scale-free, as the Plummer, the Hernquist, 
the NFW and the Gaussian profile. In all these cases we 
found that, when the initial velocity dispersion is set to zero, 
the density profiles converges, for large enough radii, at the 
same behaviors (i)-(iii) discussed above. These re sults are in 
agreement with IVisbal. Loeb fc Hernquistl l|2012h who have 
recently found that a Plummer sphere, originally at rest, 
collapses, in simulations where non-radial motions are both 
present or suppressed, forming a virialized state with a r~ 4 
tail at large radii. Note that a simple explanation fo r the 
formation of the r~ was prop osed by IWhitd (| 19871 ); Ijaffei 
{ 19871 ); in ISvlos LabmH (|2012h we introduced a more de- 
tailed model of particle dynamics, which considers as well 
the behaviour of the bound particles with almost zero total 



energy, and a series of analyses that corroborate the hypoth- 
esis that particle with very small binding energies form the 
r~ 4 tail. 

We thus conclude that when a cloud undergoes a vio- 
lent collapse , i.e. when part of its initial particles gain en- 
ergy during the collapse to eventually escape as free ones, it 
reaches a virialized quasi-equilibrium configuration that dis- 
plays properties which are formed by the dynamical relax- 
ation mechanism. In particular, these features characterise 
the tail of the density and radial velocity dispersion pro- 
files and they are related to the behaviour of those bound 
particles which were initially in the outer shells of the dis- 
tributions. 

Instead, strongly bound particles, i.e. those in the cen- 
tral dense core, have properties which depend on the initial 
conditions. These are the particles which already initially 
were strongly bound, given the cusp of the density profile: 
for this reason the energy exchange during the violent phase 
of the collapse represents a small perturbation that is not 
able to alter significantly their trajectories. 

As a final remark we note that the systems considered 
here, although pose a series of interesting theoretical ques- 
tions, differ for many important aspects not only from real 
systems but also from structures formed in cosmological N- 
body simulations. In particular, a more detailed study of the 
modification of the isolated system dynamics due to space 
expansion, continuous mass accretion, merging with macro 
objects and tidal interactions with background structures 
and a careful analys i s of n umerical and finite size issues 
Ijovce fc Svlos Labinil (|2012l ) is necessary to clarify the re- 
lation between the simple self-gravitating clouds considered 
in this paper and halos formed in cosmological simulations. 

I thank Michael Joyce for useful discussions and com- 
ments, Roberto Ammendola and Nazario Tantalo for the 
valuable assistance in the use of the Fermi supercomputer 
where the simulations have been performed. 
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